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This paper describes the advantages of having an objective classification system for 
DF bearings. The Brooke system is described in some detail, and the problems involved 
in setting up a system on these lines are considered. 



1. Introduction 

For many years we used a simple ABCD classi- 
fication system for HF DF bearings taken with 
U-A(lcock direction finders. This was not entirely 
satisfactory for various reasons. The fundamental 
difficulty was that, although in theory the classifi- 
cation depended in a well-d(4ined way on the con- 
sistency of the observed bearings, in practice it was 
largely subjective. The Brooke Variance Classifi- 
cation System is based on the Ross-Bai'field system 
developed during the Second World War and 
described by Ross [1947]. The present paper is an 
attempt to describe the system from the point of 
view of a group that might be considering the intro- 
duction of such a system. 

Section 2 of this paper reviews the reasons for 
having a bearing classification system. The sources 
of error in DF bearings are described briefly in 
section 3. The history of the Broola^ system is 
reviewed in section 4. Section 5 is (hn^oted to a 
discussion of tlu^ problems involved in iiitrochicing a 
sunilar system into another network. The solutimis 
found for the Brooke system arc indicated. Section 
6 describes the statistical analysis required to esti- 
mate variance components from check bearing data. 

2. Why Does One Need a Classification 
System? 

The first question is, wdiy does one need a beaiing 
classification system? It is impossible to give a 
short answer to this question, since it depends on 
the answer to an even more fundamental question: 
What use is one trying to make of one's DF bearings? 
One will presumably plot the bearings reported on 
any task on a gnomonic chart. One can then either 

(1) gaze hi wonder at the set of plotted lines, and 
finally emerge with a point estimate for the position 
of the transmitter, and a circle of arbitrary radius 
centered on this point such that the true position 
^'probably" lies within the circle, 
or (2) use one of the many unsound plotting methods 
that have been proposed from time to time, 

1 Contribution from Admiralty Research Laboratory, Teddington, Middlesex, 
England. 

2 Paper presented at the Conference on Transmission Problems Related to 
High- Frequency Direction Finding, at UCLA, June 21-24, 1960. 



or (3) use an electronic computer to find a more 
precise least squares solution, i.e., to find 

(a) a Best Point Estimate (or Most Probable 
Point) defined as the point minimizing the weighted 
sum of squares of the angular errors, tlie weights 
being inversely proportional to the assumed 
variances of bearing erroi's at the station concerned, 
and 

(b) a region, norniaUy a circle or a rectangle, 
approximating to the region where this weighted 
sum of squares exceeds its minimum value by less 
than some given constant, such as 4. The region 
may be called a 90 percent (say) Probability 
Region, or Confich^ice Region. Its pur[)ose is to 
give a!i a])preciation of tlu^ prol)abl(\ accui'acy of 
the point estinnite/^ 

or (4) use a method based on some valid a])])r()xiina- 
tion to the least squares solution. 

In cases (1) and (2) it may be helpful to liave some 
rough measure of the accuracy of any particular 
l)earing, but the traditional ABCD classification 
based on the consistency of repeated ol)servations, 
or some otlier sul)je('tive method, may be good 
enougli. 

In cases (3) and (4) it is important to have a 
realistic measure of the probable accuracy of any 
parliciihii- bearing. Ross [1947] wrote 

''It has long been customary for bearings to be classified 
somewhat arbitrarily by the direction-finding operator. This 
classification, intended to give an indication of the probable 
reliability of the bearings, is usually based on the quality 
(sharpness) and steadiness of the minimum. This practice 
is a legacy from the days when radio communication was 
conducted chiefly on the long and medium wavelengths, and 
it was not unreasonable to expect the operator to estimate the 
reliability of bearings in this way. In the short-wave band, 
however (20 to 100 m approximately), where propagation is 
chiefly controlled by ionospheric reflections, conditions are 
far from simple. Bearings are generally in a state of con- 
tinuous variation and it becomes increasingly difficult for an 
operator to classify them by any more or less intuitive 
process." 



3 The conventional justification for this as a Confidence Region is based on 
the assumption that all possible bearing lines from a given station are approxi- 
mately parallel. Beale [1960] provides a method of justifying it even when this 
approximation is unacceptal)le. Beale's paper shows that the approximate 
confidence region is justified if a quantity N^, called the "intrinsic nonlinearity 
of the model" is smaller than about 0.1. It seems likely that the DF pioblem 
nearly always satisfies this condition; even though one can easily imagine situ- 
ations witli a "narrow base line" when the prol)lem is decidedly nonlinear when 
exi)r('ss('<| in terms of natural parameters, such as the latitude and longitude o 
the transmitter. 
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Ross proceeds to describe an objective classifica- 
tion system which is the basis of the Brooke system. 
The details of such systems will be described in 
section 4. Meanw^iile we should consider another 
advantage in having a rational and objective 
classification system: it is almost indispensable to a 
thorough quantitative analysis of the magnitudes of 
the various errors arising in practical DF. It is 
virtually impossible to get a set of observations that 
is statistically balanced with respect to all relevant 
factors. So, for example, when studying the effect 
of distance one must allow for the difference in fre- 
quencies of transmitters at different distances; and 
such a process requires the estimation of variance 
components due to the different errors making up 
the observed bearing error. 

3. Sources of Error in HF DF Bearings 

The sources of error in HF DF bearings are dis- 
cussed by Ross [1947] and Bowen [1955]. They can 
be summarized as follows: 

(1) Instrumental errors due to small errors of 
antenna balance, and similar causes. These errors 
can be greatly reduced by a local calibration, but 
there may be an appreciable residual error varying 
with radio frequency and direction of incidence (in 
azimuth and elevation) and from site to site. 

(2) Distant-site errors, due to irregularities of, and 
obstacles on, the terrain around a high-frequency 
direction finder but outside the near site or calibra- 
tion area. The values of these errors vary with 
radio frequency and direction of incidence (in 
azimuth and elevation) . The variances vary from 
site to site, and tend to decrease as the radio fre- 
quency increases. 

From our present point of view, errors of types 
(1) and (2) are similar. They are virtually inde- 
pendent of time. For a given transmitter working 
on a given frequency the errors may vary a little 
in that the angle of elevation of the incoming ray 
will depend to some extent on ionospheric conditions, 
but there will be an appreciable systematic compo- 
nent of such an error. 

(3) Ionospheric lateral deviation errors, due to 
large scale irregularities in the ionosphere. These 
errors vary only slowly with time (with periods of 
the order of tens of minutes). Their variance does 
not depend on the DF site, but depends strongly on 
the distance of the transmitter. Not much is known 
about the effects of time of day and radiofrequency 
on lateral deviation. It seems possible that errors 
are somewhat greater at night than by day and some- 
what greater on low frequencies than on high fre- 
quencies. If so, this is probably due to the increased 
number of multiple reflections possible at night and 
on lower frequencies. Large systematic and random 
effects can be observed around sunrise. Later 
studies have not essentially changed the views ex- 
pressed by Ross [1949] on the subject of lateral 
deviation. 

(4) Wave interference errors and polarization 
errors. These errors may be quite large, but they 



generally vary rapidly (with periods of the order of 
seconds). Their effect can therefore be greatly 
reduced by taking the average bearing (preferably 
weighted according to signal strength) over a period 
of 10 sec or longer. 

4. History of the Brooke System 

The Brooke Variance Classification System arose 
out of the system developed during the Second 
World War by Mr. W. Ross and Dr. R. H. Barfield, 
and described by Ross [1947]. The basic assumption 
underlying the system is that the error in any ob- 
served DF bearing is a sum of three uncorrelated 
errors : 

(1) An instrumental and distant site error, whose 
mean square value can be determined as a function 
of the radiofrequency and the DF site. 

(2) A lateral deviation error, whose mean square 
value can be determined as a function of the distance 
from the transmitter to the DF site. This component 
must of course be estimated at the plotting center. 

(3) A wave interference and polarization error. 
An estimate of the mean square value of this error 
can be determined from the mean '^swing'' of the 
observed bearing over a period of 1 or 2 sec, or from 
the scatter of repeated observations made at intervals 
of about 5 or more sec, or both. The distance from 
the transmitter to the DF site can also be used to 
throw some light on the probable magnitude of this 
error, since wave interference due to multipath 
arrival is apt to be more serious on the more distant 
targets. 

When after the war attention was again directed 
to the classification of HF DF bearings, the Ross 
system was examined but considered to be unsatis- 
factory in its existing form. The users also con- 
sidered the subjective method of classifying bearings, 
then in use, to be unsatisfactory. Mr. Norman 
Brooke, while a member of the Admiralty Depart- 
ment of Operational Research, carried out an investi- 
gation into the problem of classifying HF DF 
bearings, the outcome of which was a proposal to 
modify the Ross system, and to extend it to apply 
to visual, as well as aural null, equipment. He 
found that it was necessary to give the variance 
components values that were rather larger tluxn 
those proposed by Ross. 

The modified version was found to be satisfactory. 
It was decided, wdth Mr. Ross's concurrence, that 
the system should be called the Brooke system. 

Before the system was introduced a statistical 
analysis was made of a large number of check bear- 
ings on known transmitters, in order to derive 
realistic estimates of all the variance components. 
Further, it is known that station performance 
changes to some extent with time, so the analysis 
of check bearing data continues on a routine basis 
in order to keep track of these changes. It is found 
that the variance component depending on the DF 
station has to be changed by 1 imit (i.e., 1° squared) 
from time to time — perhaps once every 2 years on 
the average. The other variance components are 
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reviowod IVoin time to time, but clianges have not 
been found iKM'cssaiy. 

So wo may conolude that the Brooke system works 
satisfaciorily. The essential reason for this seems 
to be that the system combines physical theory and 
empiricism in about the right pi'oportions for our 
present state of knowledge about HF DF. We 
certainly do not maintain that no further improve- 
ments in the system are possible. 

5. How Could a Similar Classification Sys- 
tem Be Introduced Into Another Organ- 
ization? 

Suppose that some authority has decided that it 
would be worth introducing an objective bearing 
classification system on these lines into liis HF DF 
network, w^iat steps must be taken to implement 
this decision? 

The first task is to define a set of factors on which 
the variance is assumed to depend. It would sc em 
natural to include 

(a) The DF station, 

(b) the radiofrequency, 



(c) the distance of 

station, and 

(d) some measure of 

bearing. 
Some discussion on 



th 



Iransiniltei* Iron 



th 



\i' quality of tlie observed 



this last item is perhaps 
required. A cursory reading of Ross [1947] miglit 
give the impression that an objective chissification 
system is only possible if one takes repeated snaj) 
observations to determine the bearing. This is not 
so. The true situation is as follows: 

Admittedly if one takes repeated snap observations 
(which may be mental averages over a period of 1 
or 2 sec), then one is likely to obtain a more accurate 
bearing than if one simply gazes at the display for 
about 30 sec, because experience has shown that it 
is almost impossible to average mentally over a 
period of more than about 5 sec. Furtliermore, the 
scatter of the individual snap bearings provides 
additional information concerning tlie reliability of 
the mean bearing. But if the transmitter is not 
active for long enough to take repeated observations, 
or if the bearing is required so urgently that there 
is not time to take and report the results of several 
snaps, then one can still use an objective classification 
system based on other factors. Indeed both the 
original Ross-Barfield system and the Brooke system 
include means of assessing a variance when it is 
impossible to obtain more than one snap. 

The next task is to postulate a formula for the 
variance as a function of these factors with a finite 
number of unknown parameters. 
It seems natural to write 

B=Vf+Va+V„ 

where B (for Brooke) denotes the total variance (i.e., 
the total mean square error if one proceeds as if all 
errors had zero means). 



Vf denotes a component depending on radiofre- 
quency and the DF station, representing instrumen- 
tal and distaiit site errors, 

Vd denotes a compouent depending on the distance 
from the transmitter, representing lateral deviation 
error, and Vq denotes a component depending on 
the quality of the observed bearing, and the type of 
equipment and display used, representing wave in- 
terference and polarization error, and also observa- 
tional error due to imperfections in the DF operator. 

We must still specify tlw forms of Vf, Vdy and Vq 
to reduce them to a finite mnnber of parameters. 
It may be helpful to consider each in turn with 
reference to the existing Brooke system. 

The Brooke system assunu^s tliat 

Vf^aj+bjf, 

where aj and bj are paranu^ers estimated separately 
for each DF station, and / denotes the frequency 
band. The definition of these frequency bands is 
being reconsidered. There is some evidence that 
the variance is not very sensitive to radiofrequency 
above about 9 Mc/s. The following bands may 
therefore prove appropriate. 

f=l for frcMiuencies between 2.000 and :^.999 Mc/s, 
f=2 for frequcMicies between 4.000 and S.999 Mc/s, 
y=3 foi' fre(|U(Micies above 9.000 Mc/s. 

For the (hstance component, the Brooke system 
assumes that 

y,= 0.6/^/2 for r/< 0.8, 
Va=l for0.8<f/<3.5, 

Va=S for:3.5<rf, 

where d ckMiotes distance from tiie transmitter to 
the DF station, measured in tliousands of kilometers. 
For the quality component, the Brooke system 
assumes that 

V,= Vs+Vn, 

where 

Vs=0 for swings between and 8°, 
= 1 for swings between 9 and 13°, 
= 2 for swings between 14 and 18°, 
=4 for swings between 19 and 23°, 
= 6 for swings between 24 and 37°, 
= 9 for swings over 38° or unmeasured, 
and Vn = r^/ndr?, 

where n is tlie number of snap bearings taken, r is 
the range of these snaps, and dn is the average range 
of n observations from a Gaussian population with 
unit variance (tabulated by Tippett [1925]). 

The theory underlying the use of these formulas is 
that Vn might be expected to give a valid estimate 
of the variance of those components of error (due to 
polarization errors and wave interference) that vary 
rapidly. But in practice the operator cannot help 
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being biased by tlie iviiowledge that all his snaps 
refer to the same transmission. He will therefore 
tend to make his observations more consistent than 
they should be, and an additional variance compo- 
nent must be introduced to allow for this. The 
bias is likely to depend on the difficulty of the task 
of taking a bearing, which is measured by the 
average ''swing". For an aural null equipment the 
swing is defined as the width of the arc about the 
minimum over which there is a just detectable change 
of signal. For any equipment with a visual display 
it is defined as the angle through which the instan- 
taneous value of the bearing moves during the time 
taken to read its mean value for an individual snap. 
There is of course no logical reason why tlie formula 
for Vs should be the same for different types of 
display, even though all the equipments have the 
same antenna system. But the statistical analysis of 
check bearings suggests that there is no important 
difference in this respect between aural null, spinning 
goniometer, and twin-channel c.r.d.f. displays.^ 

If only one snap is taken, increased values of Va 
and Vs are used to compensate for the absence of 
the Vn term.. 

Different organizations may find it convenient to 
€stim.ate V^ differently. This must depend on what 
information the DF operator finds it convenient to 
report. If he simply reports a letter classification 
(A, B, C, or D), then one could put ^^=0 for an A 
class bearing, and estimate 3 constants, being the 
values of Vq for B, C, and D class bearings 
respectively. 

One might also use the signal strength as an 
mdicator of bearing quality. This does not in fact 
enter into the existing Brooke system.. 

Having defined the form, of the bearing classifica- 
tion system., one must estim.ate the parameters. 
One can get preliminary estimates of these from 
physical theory and from, the values used in existing 
system.s, but these estim.ates m.ust be corrected (or 
confirm.ed) using data obtained by the DF organiza- 
tion in question. These data m.ay be either genuine 
operational data or check bearing data where the 
true positions of the transmitters are known to the 
analyst. 

As far as I know, no serious attem.pt has yet been 
m.ade to estim.ate variances from, operational data 
when the true positions of the transmitters are 
unknown. The use of such data poses the funda- 



< There is a small point of some statistical interest about the estimation of Vq. 
Elenientary theory suggests that V Vn is an unbiased estimator of the standard 
deviation of errors varying randomly from snap to snap, for a given value of this 
standard deviation. But this does not prove that V Vn provides the best estimate 
of t his standard deviation, given the range r. Indeed one might expect that a 
belter estimate would be obtained by attaching some weight to the estimate 
obtainiMJ iioiii the long-term average value of this standard deviation— thereby 
makin'j I he used estimate less dependent on r. After holding forth on these lines 
to niv cMiJcaiiueson more than one occasion, I eventually persuaded them to do a 
delailed analysis of the check bearing data to investigate this point. They 
found that, if anything, the variance component Vn should be made more de- 
pendent on r. It has been suggested tliat this is because the dividing factor 
dn'^ is too large because the error distiibution is not normal. (And indeed even 
from a normal population a sniallei- dividing factor is required to give an unbiased 
estimate of the variance as opposed to t he standard deviation.) But I prefer to 
regard this as suggesting that the iang(> has additional value as an indicator of the 
quality of the bearing. More light could be thrown on this point by an analysis 
of the "individual snap observations, but it is of little operational importance. 



mental difficulty that one has to lean rather heavily 
on the assum.ption that the bearing errors are uncor- 
related. If one is prepared to accept this assump- 
tion, then it is fairly easy to estim.ate the absolute 
values of the variances of a set of DF stations if we 
know their relative variances. For in this case we 
can determine the true Best Point Estimate; i.e., the 
hypothetical position of the transmitter m.inimizing 
the sum. of squares of the bearing errors divided by 
their variances. We can then use the fact that (if 
the error distributions are norm.al with zero m.ean) 
the m.inim.um. value of this sum. is distributed as x^ 
with n-2 degrees of freedom., where n is the num.ber 
of DF stations contributing to a task. In particular, 
the mean value of this sum. is n-2.'^ 

In practice the relative variances are unlikely to 
be known. The problem, then becomes considerably 
more complicated, ft is discussed in a companion 
paper, Beale [1961]. 

It is much easier to estim.ate variances from, check 
bearings, provided satisfactory check bearings can be 
obtained. 

Check bearings m.ay be used both to verify that 
the equipment is functioning; properly and to indicate 
the probable accuracy of the DF station on opera- 
tional tasks. A check bearing program, to check the 
equipment does not need to be controlled nearly as 
carefully as a clieck bearing program, to estim.ate 
accuracy. Bo wen [1955] writes 

''The use of check bearings on known transmitters in 
establishing the performance of a high-frequency direction- 
tinder is a widely accepted technique. It has been found, 
however, that two problems present themselves: 

(a) There is difficulty in arranging a programme in which 
the radio frequencies, ranges, azimuths, propagation paths, 
transmitter powers and other parameters are adequately 
sampled and distributed similarly to those w^hich will define 
the normal task of the direction finder. 

(b) Unless great care is taken to hide the identity of the 
transmission from the operator, both as regards its true 
bearing and its identity with previous transmissions, con- 
siderable operator bias will occur. 

"The methods by which these problems are overcome will 
generally be peculiar to the task with which the d.f. organi- 
zation is concerned. It may be that (a) will require the 
deliberate setting-up of a transmitting system, although care- 
ful selection of known trans .nitters may obviate this; (b) can 
be overcome either by careful organization of the way in 
which check-bearing tasks are fed to the operator, or by a 
system of random scale displacements at the direction-finder.'' 

Even if one fails to overcome these problems 
completely, it would seem unreasonable to use this 
as an excuse for not introducing an objective classi- 
fication system. Indeed the results from an objec- 
tive system should be very much better than those 
derived from overall station performances based 
on such data — since the system will give some 
rational method of extrapolating to conditions met 
only rarely in the check bearing program. 

The statistical problems involved in estimating 
variance components from check bearings are dis- 
cussed in the next section. 

£ Provided, as will usually be the case, that the quantity N<^ introduced by 
Beale [1960] is small. 
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Estimating Variance Components 
Check Bearing Data 



From 



The ruiulaiiuMital statistical problem in this work 
can be expressed as follows: 

Given // independent observations, each from a 
normal population with zero mean, such that the 
it\i observation Xi has a variance 



Br- 






(6.1) 



where the atj are all known but some at least of the 
Vj are unknown, estimate the unknown Vj. 

Before showing how to solve this problem, I 
should perhaps indicate its relevance. The quantity 
Xi represents the actual error in the -ith bearing 
observation, Bi represents its assumed variance, and 
the Vj represent the parameters in the classification 
system. For example, suppose that we have a 
simple model with one component that can be 
assumed known, and unknown (^oiuponents depend- 
ing on the radio frequency and quality, assumed to 
depend on 4 parameters, such that 

Vf=Vi+Jv2, 

where y denotes frequency band, eitlier 1, 2 or 3, and 

Vg = Vs if the swing <8°, 

Vq—V4, if the swing >8°. 

Now suppose that the first few check l)earings 
have the following characteristics: 

Known Vari- Frequency 
Bearing N'o, ance Coinjponent Band Swing 



1 


6 


1 


3°, 


2 


3 


2 


13°, 


3 


5 


3 


7°, 



etc. 

Then, if we write formally 2^5= 1, (6.1) is satisfied 
with 

aii = l, 0:12=1, 0^13^1, 0^14 = 0, «15 = 6, 
a2i = lj Oi22 = 2, ^23 = 0, 0:24=1, <^25 = 3, 



0^31 



= lj 0:32 = 3, 033=1, 034 = 0, 035 = 5, 



the parameters to be estimated being Vi, V2, v^^ and ^4. 

We assume that our bearing errors are measured 
ni degrees. Variances are therefore measured in 
degrees squared. 

Now 7?i is intended to represent the mean value 
of xi, so it seems natural to try to minimize the sum 
of squares of the deviations of the xi from the corre- 
sponding Bi. But we must remember that the ex- 
pected magnitudes of these deviations depend on 
the variances of the xt In fact if Xf is normally 
distributed the variance of a^f is 25|. It therefore 



seems reasonable to try to minimize 



i=l 



(6.2) 



where the weights Wf are regarded as constants but 
where their numerical values are given by 



Wi=l/Bl 



(6.3) 



Differentiating (6.2) witli respect to Vj, regarding 
Wi as constant, we have, using (6.1) 



^Wiaij{x1 — Bi)=0, 



(6.4) 



for each^* such that Vj is to be estimated. 

In particular if for some j, a^ is always either 
or 1, the corresponding equation from (6.4) 
becomes 

S>,(x?-5,) = 0, 

where 2- denotes summation over all observations 
such that 0:^^=1. 

The problem is therefore to find estimates Vj for 
the unknown Vj satisfying (6.1), (6.3) and (6.4). 
This will have to be done iteratively. Choose trial 
estimates for the Vj, compute the corresponding 
values of Bi from (6.1), and hence compute the 
weights Wi from (6.3), substitute these values of 
Wi in (6.4) and hence estimate new values of the v,.^ 
The whole cycle can be repeated if necessary, but 
it should be remembered that snudl (relative) errors 
in the Wi will be of little significance. 

It can easily be shown that these estimates are in 
fact maximum likelihood estimates for the Vj. The 
above derivation is less rigorous but seems more 
intuitive.^ 

One important proviso should be made. It is 
undesirable to have weights for individual bearings 
of the same order of magnitude as the sum of the 
weights for all otlier bearings put together. In the 
analyses for the Brooke system, this is avoided by 
putting Wi=y2 whenever Bi<Cr<J2, So the used 
version of (6.3) reads 

^,=min (K, IIB\), 

It is of interest to consider the probable ac- 
curacy of the final variance estimates. It is easy 
to find the variance of an individual component, 
sa}^ Vk, if we assume that all the other components are 
known exactly, and also that each Wi equals the 
reciprocal of the square of the true variance of the 



« In practice it may be more convenient to improve the estimates of the param- 
eters referring to one variance component, say Vf, assuming that the others 
have their present trial values. This reduces the number of equations that have 
to be solved simultaneously. 

7 Some statisticians may point out that for this problem the method of maxi- 
mum likelihood gives an "inadmissible" estimator, since it may produce negative 
variance components. With a fair-sized sample this is unlikely to happen; but, 
if it does, one must replace the negative component by some more plausible value 
and re-estimate the other components accordingly. 
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corresponding Xi. For (6.4) then becomes 
so that the estnnate Vjc of Vjc is given by 



i=l ;?^/c 

Now var x'i=2Bj^ and Wi=\IB\, so 



var t^;t 



^±(||J..«/(S|J 



=2/Sw,aL 



If Qf,t=0 or 1 for all i, (6.5) reduces to 



var ft=2/2>(, 



(6.5) 



(6.6) 



where S;^ denotes summation over those check 
bearings for which a^A;— 1, i.e., for which the given 
variance component is relevant. 

Errors in estimating the weights Wi^ and also 
errors in the other variance components, will not 
normally have a major effect on the variance of 
Vie* On the other hand the fact that the observa- 
tions may not be strictly independent, and that the 
error distribution may not be strictly normal, may 
have a considerable effect on the validity of these 
formulas. Our experience suggests that, presumably 
because of these causes, formulas (6.5) and (6.6) 
underestimate the variances of the f^A; by a factor of 
about 2. The formulas 

var Vk^^l^Wia^k, or var Vk-^l^[Wi, 

are therefore recommended in practice. 

There is one pecularity about these variance 
component estimates that one should be aware of, 
although it does not affect the operation of the system 
in any way. There is a degree of indeterminacy 
about the system in that a constant can be sub- 
tracted from all values of, say, Va and added to all 
values of, say, Vf without altering any of the Brooke 
variances. The only restriction is that none of the 
variance components can ever be negative. Our 
practice is to put the minimum values of Va and Vq 
equal to zero, and to put the '^unattached variance" 
into Vf. This is purely for operational convenience. 
In fact we believe that the minimum value of the 
effect of lateral deviation, which is represented by 
Vd, is about 1° squared, as indicated in section 5 
above. 

We have now considered the estimation of variance 
components in the first instance. 



The other vital element in an objective bearing 
classification system is a method of keeping the 
system up-to-date. For routine corrections to be 
timely and not unduly subject to sampling fluctua- 
tions, the corrections should be based on as simple a 
statistical model as possible. We therefore assume 
that all variances of bearings taken at a particular 
station are underestimated by a constant small 
amount. The best estimate of this amount is then 
approximately 

^Wi {x',-Bi)/i:Wi, 

where summation extends over all bearings taken 
at this station. This quantity is called the ' 'apparent 
correction'^, short for ''apparent Brooke variance 
correction", audits variance is approximately 4:jl^Wi. 
If we wait until S'K;f>20, the variance of the 
apparent correction is less than about 0.2, so the 
correction is unlikely to be in error by as much as 
one unit. Nevertheless the apparent correction is 
rounded down in absolute value to the nearest unit, 
since some weight should be given to the fact that^ — ■ 
if the system has been running for some bime^ — 
previous data suggested that no correction was 
necessary. 

In practice this works as follows. For each check 
bearing the ^'weight" Wi and the "indicator" 
'^i{xi — Bi) are recorded, and these are summed for 
all check bearings taken in a month at each station. 
These data are inspected and combined with those 
for enough previous months to make l^Wi> 20 for the 
station in question. Then l^Wi(x^i—Bi)l2Wi, rounded 
down in absolute value to the nearest unit, is added 
to all values of Vf for this station; i.e., to the value 
of ttj in the formula Vf=aj-\-b{f. 

After any adjustment to the value of aj, the 
monthly totals of l^Wi{x\—B^ must obviously be 
reduced by 1>Wi multiplied by the addition to a^ 
before they are combined with subsequent monthly 
totals to decide whether further changes should be 
made. But it should never be necessary to re- 
compute the Bi and Wi for individual bearings. 

In practice our estimated variance components are 
not subject to frequent oscillations. 

In addition to these routine adjustments to the 
overall variances for each station, it is important 
that the other variance components should be re- 
viewed from time to time, by computing apparent 
corrections from the formula 'EWi(x'i—Bi)/XWi with 
summation extending over all check bearings asso- 
ciated with a particular value of the component 
under review. The desirability^ of introducing fur- 
ther terms into the expression for the variance can 
be tested in the same way. 

This may be as good a place as any to mention 
the problems of systematic errors and wild bearings. 
Many DF stations exhibit small systematic errors 
in their bearings. These have been studied for some 
time, but corrections to allow for them are not ap- 
plied because the physical basis of the errors is not 
understood. If the systematic error is ignored, one 
is acting as if the error distribution had mean zero. 
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and the Brooke variance is really an estimate of the 
mean sqnarc error. 

DF stations are apt to produce a certain number 
of wild l)earings that are a long way off the true 
bearing, and which one hopes to be able to reject 
in plotting by their inconsistenc}^ with the other 
bearings. Obviously one does not want to include 
these in one's check bearing data for the determina- 
tion of variances. But the problem is, how large 
must the error be before one is justified in rejecting 
the bearing as wild on the grounds that it would 
probably be recognized as such even when the true 
target position was unknown? This is not an easy 
problem. The proper answer presumably depends 
on the variance of nonwild observations, and on the 
number, location, and accuracy of the other DF 
stations in the organization working on the same 
tasks. The rejection of bearing errors greater than 
10 or 15° may often be reasonable. 



Although it has fallen to me to present this account 
of the Brooke system, I liope I have made it clear 
that jnost of the crcnlit l)elont!:s elsewhcM'e. I have 



been abh^ to refer to some of the relevant work by 
Messrs. Ross, Brooke, and Bowen; but many people 
have contributed in various ways — ^notably tliose 
who organize the check bearing program, and collate 
and analyze the data. 

7. References 

Beale, E. M. L., Coiifidonce regions in nonlinear estimation, 

J. Roy. Stat. Soc. (B) 22, 41-88 (1960). 
Beale, E. M. Ij., The estimation of variances of position 

lines from fix(\s with unknown target positions, J. Research 

NBS 65D (Radio Prop.) (May-June IDtU). 

Bowen, K. C, Sources of error in U-Adcock hi^h-freqiiencv 

direction-finding, Proc. I.E.E. 102 Pt. B, 529-532 (1955). 
Ross, W., The estimation of the probat^Ie accuracy of high 

frequency radio direction-finding bearings, Journal I.E.E. 

94 Pt. 3 A, 722-726 (1947). 
Ross, W., Lateral deviation of radio waves reflected at the 

ionosphere, DSIR Radio Research Special Rept. No. 

19 (1949). 
Tippett, E. H. C, On the extreme individuals and the range 

of samples taken from a normal popidation, Jiiometika 17 

364-387 (1925). 

(Paper 65D3^128) 



261 



